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ABSTRACT 

The canonical statistical model analysis of strange and multistrange hadron production in 
central A-A relative to p-p/p-A collisions is presented over the energy range from yfs = 8.73 
GeV up to \/s = 130 GeV. It is shown that the relative enhancement of strange particle yields 
from p-p/p-A to A-A collisions substantially increases with decreasing collision energy. It is 
largest at y/s = 8.7 GeV, where the enhancement of Q, S and A is of the order of 100, 20 and 
3, respectively. In terms of the model these results are due to the canonical suppression of 
particle thermal phase space at lower energies, which increases with the strangeness content of 
the particle and with decreasing size of the collision fireball. The comparison of the model with 
existing data on energy dependence of the kaon/pion ratio is also discussed. 
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1 Introduction 



The production of strange particles is extensively studied in heavy ion experiments in a very 
broad energy range from GSI/SIS to BNL/AGS and from CERN/SPS up to BNL/RHIC 0. 
The amount of data available for K + and K~ yields already allows a detailed analysis of the kaon 
excitation function in A-A collisions [Q, |3|, || as well as an analysis of the relative enhancement 
of kaon production in A-A with respect to p-p collisions || || || . Of particular interest is here 
the behaviour of K + yield with energy, as this is the most abundantly produced particle with 
non-vanishing strangeness. The kaon excitation function at midrapidity, as the K + /n + ratio, 
is a very abruptly increasing function of the collision energy between SIS up to top AGS. At 
higher energies it reaches a broad maximum between 10 AGeV and 40 AGeV || and, gradual 
decrease up to RHIC energy |7], In the microscopic transport models the increase of the 
kaon yield with collision energy is qualitatively expected, owing to a change in the production 
mechanism from associated to direct kaon emission. However, the hadronic cascade transport 
models are as yet not providing the quantitative explanation of the experimental data. The 
Hadron String Dynamics (HSD) model |§ severely underpredicts the top AGS results on K + /tt + 
ratio, while RQMD || overpredicts the yield at lower energies and gives too small yield at the 
SPS. The Statistical model SM, on the other hand, provides quite a satisfactory description 
of the kaon excitation function and of the K/ir midrapidity ratio in the whole energy range 
from SIS up to RHIC (T^, [□]]. A detailed analysis of the experimental data at SIS jnj, AGS 



1^, |l3fl , SPS [|T3|, [L4| and recently also at RHIC |0| [16], |I7| has shown that hadronic yields and 
their ratios resemble those of a population in chemical equilibrium. Most particle multiplicities 
measured in nucleus-nucleus collisions are well consistent with the thermal model predictions 



13, 14, 15, 18]. A similar analysis of particle yields in hadron-hadron collisions [19] also shows 



the consistency of the statistical model with the data. 

Recently, it was argued p0[ that SM can also describes a basic future of strangeness en- 



hancement from p-A to A-A collisions measured by WA97 Collaboration [21]: the enhancement 



pattern and enhancement saturation for a large number of participants A^ part . The enhancement 
of (multi) strange baryons is of particular interest, since it was conjectured to be a signal for 



quark-gluon plasma formation in heavy ion collisions [2~2"|. The centrality dependence of the 
(multi) strange baryon enhancement, increasing with the strangeness content of the particle, 
was shown to appear in the SM as a consequence of the canonical suppression of the particle 
phase space in a small system. Satisfactory agreement of the SM with data confirmed that SM 
provides a tool to make a prediction for hadron production in A-A collisions. 

In this paper we apply (SM) formulated in the canonical ensemble [^, |24| to study the 
energy dependence of the (multi) strange baryon enhancement from p-p to A-A collisions. We 
first show that the observed relative enhancement of the K + /tt + ratio in A-A to p-p collisions 
from top AGS up to top SPS energy can be quantitatively described by the SM. We then 
make a prediction for relative (multi) strange baryon enhancement at ^/s = 8.73, 12.3, 17.3 GeV 
at SPS and ^Js = 130 GeV at RHIC energies. We show that at RHIC the enhancement of 
(multi) strange baryons is expected to be comparable with the one measured at the top SPS 
energy; however, it could be almost an order of magnitude larger at lowest (y^ = 8.73 GeV) 
SPS energy. 



1 



2 Statistical description of strangeness enhancement 



Within the statistical approach, hadron production is commonly described by using the grand 
canonical GC partition function, where the charge conservation is controlled by the related 
chemical potential. In this description a net value of a given U(l) charge is conserved on 
the average. The GC approach can be valid only if the total number of particles carrying 
a quantum number related with this symmetry is very large. In the opposite limit of small 
particle multiplicities, conservation laws must be implemented exactly and locally, i.e. the 
canonical C ensemble for conservation laws must be used |23|, |26|. The local conservation of 
quantum numbers in the canonical approach severely reduces the thermal phase space available 
for particle production. The exact charge conservation is therefore of crucial importance in the 
description of particle yields in proton-induced processes |TP|, [2(|, in e + -e~ as well as in 
peripheral heavy ion collisions 



In the present analysis of the energy dependence of the strangeness enhancement, we adopt 
the statistical model, which has previously been used to analyse hadron productions in heavy 
ion collisions [^, ^4j. This model is formulated in the canonical ensemble with respect to 
strangeness conservation, whereas baryon number conservation is controlled by the chemical 
potential. Within this framework particle multiplicity ratios depend only on two thermal pa- 
rameters: the temperature T and the baryon chemical potential fis, as well as the correlation 
volume, which is assumed to scale with the number of participating projectile nucleons. 

The canonical model description of strangeness conservation including multistrange particle 
contributions was described in detail in |20], [24], 27]. We quote here the final results 0] for 



particle density, which are relevant to our discussion of the energy dependence of the strangeness 
enhancement . 

The number density nf of particle i carrying strangeness s is found to be 



oo oo 



where 



n\ = —±- ala n 2 a x 2n 3p s I n (x 2 )Ip(x 3 )I- 2 n-3p-s(xi), (1) 

5=0 Ti= — oo p— — oo 

oo oo 

Z% =0 = E J2 ala2a^ 2n ' 3p I n (x 2 )I p (x 3 )I^2n-3p(x 1 ) (2) 

n=— oo p=— oo 

is the canonical partition function for a system with total stangeness 5 = 0. The variables 
a, and Xi are defined by dj = ^Si/S-i, Xi = 2^SiS~i, and S n = VJ2k %k is the sum over all 
particles and resonances carrying strangeness n. For a particle of mass rrik, with degeneracy 
factor gk, carrying baryon number Bk with the corresponding chemical potentials fiB, the one- 
particle partition function in Boltzmann approximation reads 

Zi = ^mlTK 2 (^) exp(^ B ); (3) 

here Ij and K 2 are the modified Bessel functions. 

It can be shown that, in the limit of large x\, that is for large strangeness multiplicity, the 
above formula coincides with the grand canonical result In the opposite limit, that of small 

1 This can be easily seen for charges and baryons neutral system by making a large argument expansion of 
the Bessel functions in Eq. (nl). 
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xi, the equilibrium density of (multi) strange baryons is strongly suppressed with respect to 
their grand canonical value [[2(J . 

To study the energy and centrality dependence of (multi) strange particle yields in terms 
of the above model, one needs to establishes first the variation of thermal parameters with 
energy and centrality. From a previous analysis |28] one knows that temperature to a good 
approximation is only a function of the collision energy and is independent from the number 
of participating nucleons. For fully integrated, 47r particle yields, one can also assume that the 
baryon chemical potential is weakly changing with centrality. The above assumption could be 
partly justified by data since for ^/s = 8.73 GeV the total number of pions per participant 
in central Pb-Pb ((ir) / (N pait )=2.5) and in p-p ((tt) / (N paxt ) ~ 2.62) almost coincide (see 
also Fig. 1). In the SM and with the same temperature, the above can only be valid if the 
baryon chemical potential is weakly centrality-dependent. For larger collision energies, such 
as at RHIC, (ir) / (N part ) is seen in Fig. 1 to increase from p-p to central Au-Au collisions. 
Thus, here \xb is decreasing with increasing centrality. With a freeze-out temperature T ~ 175 



MeV and //# ~ 50 MeV required to reproduce the RHIC data |2||, this would correspond to 
an increase of \ib from 50 MeV at most central Au-Au collisions to 70 MeV in p-p collisions. 
This change, however, will only weakly influence the relative enhancement of (multi) strange 
particles. We thus assume, for simplicity, that the baryon chemical potential is centrality- 
independent.0 Thermal parameters are, however, very sensitive to collision energy. At the AGS 
|T^ 1 and SPS [14| we use the chemical freeze-out parameters as obtained from a detailed analysis 
of the experimental data. In particular at a/s = 17.3 GeV, T ~ 168 MeV and fis ~ 266 MeV. 
At RHIC we adopt the values from At 40 AGeV the available data are still not sufficient 
to make a precise determination of freeze-out parameters [3C]. To estimate these values we 
use the relation between T and [ib as determined by the unified freeze-out conditions of fixed 
energy per particle ]32| and the measured result on (7r)/V r part = 2.55 ±0.15 for Pb-Pb collisions 
at 40 AGeV g. The Pb-Pb data at 80 AGeV for (7r)/A part are still not known. Here we use 
instead the extrapolated to 80 AGeV value from Fig. 1, leading to (7T)/N pait ~ 3.4. Thermal 
parameters for 40 and 80 AGeV extracted in this way are: T = 145 ± 5 MeV, ^b — 370 ± 20 
MeV and T ~ 152 MeV, ^b ~ 280 MeV correspondingly. Finally, the volume parameter 
appearing in Eq. ([[]) is assumed to be proportional to the number of projectile participants, 
V ~ Vo Ap art /2 where Vo — 7 fm 3 is taken as the volume of the nucleon. 

In Fig. we show the compilation of the data on the K + /ir + ratio in A-A relative to p-p 
collisions |3], |J. This double ratio could be referred to as strangeness enhancement factor. The 
enhancement is seen in the data to be the largest at the smallest beam energy and is decreasing 
towards higher energy. The line is a smooth interpolation between the canonical model results 
for y/s = 17.3, 12.3,8.73,5.56 GeV, obtained with the parameters described above. 

The enhancement seen in Fig. |2| is entirely due to the suppression of the K + /tt + ratio in p-p 
collisions with decreasing energy and not due to a dilution of this ratio by excess pions in the 
A-A system. The K + /tt + ratio is known experimentally not to vary within 30% in the energy 
range from y'i ~ 5 GeV at AGS up to y/s = 130 GeV at RHIC P, II, iT|. This behaviour is 



also described by the statistical model as a consequence of particular distributions of thermal 



2 In general in p-p collisions the baryon number should also be treated canonically [O. However, since here 
the baryon number B — 2 the canonical suppression is much less effective than in the case of strangeness. The 
canonical corrections due to baryon number should not exceed 30% [BjJ. 
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Figure 1: Total number of pions per partici- 
pant ((tt) / N paxt ) , 4n ratio versus energy. The 
data at lower energies in A- A as well as in p- 
p collisions are from [2{J. The RHIC results 
are from |33| . The short-dashed and dashed 
lines represent the fit to the data. 



Figure 2: K + /tt + ratio in A+A relative 
to p+p collisions. For the compilation of 
data, see 0, |[ [Xlj . The dashed line rep- 
resents the statistical model results. 



parameters with collision energy [24].f] 

The results in Fig. show that in terms of the SM the energy dependence of the strangeness 
enhancement could be related with the suppression of the thermal phase space available for 
strangeness production in p-p collisions. Having in mind the above agreement, we apply this 
model to analyse the enhancement for multistrange baryons. Two relevant questions are here of 
particular interest: (i) Is the enhancement pattern observed by WA97 a characteristic feature 
of the top SPS energy or can it also be seen at lower energies ? (ii) Is the enhancement for 
multistrange baryons a decreasing or increasing function of energy? In Fig. ^ and Fig. f| we 
show the results on relative (multi) strange baryon enhancement from p-p to Pb-Pb collisions 
at -y/i = 8.73 GeV and at = 130 GeV, respectively. It is clear that the same enhancement 
pattern as at the SPS is expected in the SM model to appear for all relevant energies. 

To see the dependence of the strength of the enhancement with energy, we show in Fig. 
|5| and Fig. |6] the relative enhancement of H and (Q + Q) for different collision energies. The 
enhancement is the largest at lowest energy; for Q, it can even be larger by a factor of almost 
10 at 40 A GeV than observed at the SPS. This behaviour could already be partly deduced 
from the experimental data. Indeed, combing the Si-Pb results for S~ production obtained by 
the E810 Collaboration |38| and the E802 value for pion or K~ yield in Si-Au collisions |39| at 



3 We have to point out that the NA49 data in the full phase space shows a depletion in the K + /ir + ratio 
from 40 A GeV to 160 A GeV. This behaviour was previously conjectured as a possible signal of deconfincment 
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Figure 3: Centrality dependence of 
the relative enhancement of particle 
yields/participant in central Pb-Pb to p- 
p collisions at fixed energy y/s = 8.73 
GeV. 



Figure 4: Centrality dependence of 
the relative enhancement of particle 
yields/participant in central Pb-Pb to p- 
p collisions at fixed energy ^fs = 130 
GeV. 



the top AGS energy, one can estimate that H /ti + ~ 0.0076; this within errors, coincides with 
the value of E-/n + ~ 0.0074 obtained by NA49 in Pb-Pb at y/s = 17.3 GeV g|, as seen in 
Fig. [F[ In p-p collisions the 5~/7r + ratio is obviously a strongly decreasing function of energy. 
From the above one could therefore expect that the relative enhancement of S~ from p-p, p-A 
to A-A collisions should be larger at AGS than at SPS energy, which is seen in Fig. |5|. In 
addition the ratio containing only newly produced particles such as e.g. H~ /K~ , is also seen in 
Fig. ^ to be larger at AGS than SPS. The energy dependence of the E~ /n + and E~ /K~ ratios 
was calculated in Fig. [7] along the freeze-out curve |32] following the method described in |24 



The results of Fig. |5] and Fig. |5] show a saturation of the enhancement, which indicates 
that the grand canonical limit was reached. It is clear from these figures that this saturation 
is shifted towards larger centrality with decreasing energy. 

In p-p collisions local strangeness conservation required that there be associated production 
of particles, e.g. 5 has to appear together with two other strangeness 1 particles, to neutralize 
strangeness. In high energy central A-A collisions, there are already sufficiently many strange 
hadrons being produced for strangeness to be conserved on the average, which substantially 
increases thermal particle production. The associated production and locality of strangeness 
conservation is an origin of the suppression of particle thermal phase space, which increases 
with the strangeness content of the particle as well as with decreasing collision energy. As a 
consequence the enhancement from p-A to A-A collisions with decreasing energy is seen in Fig. 
H and Fig. ^ to be stronger for Q than for S. At RHIC the freeze-out temperature was found 
to be consistent, within errors, with the one at the SPS as well as with the critical temperature 
found in the Lattice Gauge Theory. As a consequence the enhancement of (multi) strange 
baryons at y/s = 130 GeV RHIC energy is seen in Fig. 4 to be comparable with that at 
the SPS. The change in baryon chemical potential from 266 MeV to 50 MeV, evidently does 
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Figure 5: Centrality dependence of relative 
enhancement of S~ yields/participant in cen- 
tral Pb-Pb to p-A reactions at different col- 
lision energies. 



Figure 6: Centrality dependence of relative 
enhancement of (f2 + f2) yields/participant 
in central Pb-Pb to p-p reactions at dif- 
ferent collision energies. 



not influence the strength of the enhancement. Thus, also an increase of the energy from 
130 to 200 GeV should not change the above results much. These predictions are in obvious 
disagreement with UrQMD [|37|]. There, the production of strangeness is very sensitive to the 
initial conditions. In UrQMD the early stage multiple scattering may imply an increase of the 
colour electric field strength. Consequently, according to the Schwinger mechanism, this should 
increase the production of (multi) strange baryons. Under similar kinematical conditions as at 
the SPS, the UrQMD model predicts at RHIC an increase in relative strength of Q yields from 
p-A to A-A by a factor of 4 |37| . 

The strangeness enhancement pattern measured by the WA97 Collaboration at the SPS 
21]] was predicted as a signal for quark-gluon plasma formation p2| |. In the context of the 



considered SM model the enhancement pattern of (multi) strange baryons should be observed 
at all SPS energies, with increasing strength towards lower beam energy. Thus, the results of 
the above SM make it clear that strangeness enhancement and enhancement pattern are not a 
unique signal of deconfinement, since these features are expected to be there also at energies 
where the initial conditions are very unlikely for deconfinement. 

The quantitative results shown in Figs. contain some uncertainties. The magnitude 
of the enhancement is very sensitive to the temperature taken at a given collision energy. 
Changing T by 5 MeV, a typical error on T in a thermal analysis, can change, for instance, 
the enhancement of Q shown in Fig. |3| by a factor of 2. The N pauTt dependence of the strange 
hadron enhancement seen in Figs. 3-6 was obtained by assuming a linear dependence of a 
volume parameter V in Eq. with A" part . In general, V could have a weaker dependence 
with centrality, which could be reflected with an only moderate increase of the enhancement, 
with centrality and saturation appearing at larger volume. In addition, including the variation 
of the thermal parameters, in particular of the baryon-chemical potential, with centrality, or 
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Figure 7: Ratios of E~/tt + and S~/K~ as a function of collision energy. The data at the 
SPS are fully integrated NA49 results. The corresponding ratio at the top AGS was obtained 
from E810 results on S~ measured in Si-Pb collisions in the rapidity interval 1.4lyl2.9 |38|] 
normalized to the full phase space values of 7r + and K~ yield obtained in Si-Au collisions by 



E802 §5] 



extending the model to a canonical description of baryon number conservation [|T^], or finally 
including a possible asymmetry between the strangeness under saturation factors in p-p and 
A- A collisions 



13, 19 



could change our numerical values. However, independently of these 
the main results: (i) enhancement decreasing with increasing collision energy, 



uncertainties 

and (ii) enhancement pattern being preserved at all SPS energies, are always valid. 



3 Summary and conclusions 

In conclusion, we have shown that, in terms of the statistical model, the relative enhancement of 
(multi) strange baryons from proton-proton or proton-nucleus to nucleus-nucleus collisions is a 
decreasing function of the collision energy. Experimentally this fact already, had been obtained 
for kaon yields and is shown to be expected for multistrange baryons. In addition, an increase of 
the enhancement with the strangeness content of the particle is a generic feature of our model, 
independent of the collision energy. On the qualitative level the only input being required 
in the model to make the above predictions is the information that freeze-out temperature is 
increasing with collision energy and that the chemical potential shows the opposite dependence. 
The above conditions are well confirmed by a very detailed analysis of particle production at 
different collision energies. 

We have presented the quantitative predictions for a relative enhancement of A, 5 and Q 
yields in the energy range from y/s = 8.7 up to y/s = 130 GeV. We have discussed the possible 
uncertainties of the presented results. The relative enhancement at RHIC was found to be 
lower than at the SPS. This is in contrast with the UrQMD predictions, which are showing the 
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opposite behaviour |37j • The statistical model applied here is also very unlikely to be capable 
of explaining an abrupt change of the S/iVp art enhancement recently reported by the NA57 
Collaboration [[|3[ . 
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